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ABSTRACT 

We study relativistically expanding electromagnetic fields of cylindrical geometry. The 
fields emerge from the side surface of a cylinder and are invariant under translations 
parallel to the axis of the cylinder. The expansion velocity is in the radial direction 
and is parametrized by v = R/{ct). We consider force-free magnetic fields by setting 
the total force the electromagnetic field exerts on the charges and the currents equal to 
zero. Analytical and semi-analytical separable solutions are found for the relativistic 
problem. In the non-relativistic limit the mathematical form of the equations is similar 
to equations that have already been studied in static systems of the same geometry. 

Key words: gamma-rays: bursts; methods: analytical; MHD; stars: magnetic fields; 
stars: neutron; 



1 INTRODUCTION 

The presence of magnetic fields in objects of astrophysical 
interest is evident. Starting from the Earth's magnetic 
field, to the solar magnetically driven phenomena and 
at the most extreme from magnetars to the recently 
observed filaments in the Perseus cluster of galaxies 
l| Fabian et al.l |2008| ). magnetic fields play an important 
role in the evolution of the object they occur in. The 
study of the magnetic fields is complicated due to the 
form of partial differential equations that have to be 
solved and their non-linear behaviour. For this reason 
most of the research has been done on equilibrium con- 
figurations. Even in the non-relativistic li mit a few of the 
publis h ed studies are ana l ytica, l, for instance iPriest fc Milne 
1980| ) : lTsinganosl (I1981|): IlowI iI 19821) :lLvnden-Bell fc Boilv 



19941 ): lAlvl (l994l ): IVlahakis fc Tsinganoj (Il998l) whereas 



most o f them are c o mputa ti onal i.e . Wolfson fc Lowl 
l|l992h : iLovelace et all l|2002l ): iFendt I (|2008h and ref- 
erences therein. In the relativistic case the studies 
are mainly motivated by pulsar magnetospheres, 7- 
ray bursts, and relativistic j ets. Most of the so lution s 
are numerical, for instance [LcBlanc fc Wilson' ('l970l); 

Konigl 



IContopoulos. Kazanas fc Fendtl ( 19991: Vlahakis fc 



2004): iGoodwin. Mestel. Mestel fc Wright] 



Komissarov 



Uzdenskv fc MacFadvenI \2m% 



Meliani et al.l 



2004 ) 



20061) 



McKinney fc NaravanI 



I 2007t ): I Kalapotharakos fc C ontopoulos and only 

a few of them are analytic Lpvojacc ( 1976 ) ; iPrendergastI 
l|2005h : lGourgouliatos fc Lvnden-Belll (I2OO8I ). 

In this paper we study the problem of a force-free rel- 
ativistic electromagnetic field inspired by the studies previ- 



ously made by lAlvl l| 19941) for the geometry and the sym- 
metry used; and from IPrenderg ast (2 0051) fo r the relativistic 
generalisation we make. ^Pr endergastl ([2005') showed that it 
is possible to study magnetic fields that expand uniformly 
within a sphere. These fields emerge from a point, obey the 
ideal MHD relation E = —v x B and the net electromag- 
netic force is zero, then he solved the partial differential 
equation he found assuming a linear form. In this paper we 
show that an analogue case that permits analytical solutions 
and allows us to include time evolution exists for cylindrical 
magnetic fields. The field emerges from a linear singularity 
as in Alv, (1994) who studied this type of fields in the static 
limit, but in our problem the field expands radially. The 
expansion velocity is proportional to the distance from the 
axis and inversely proportional to time, thus every element 
of the plasma moves at constant speed. The expansion is 
slow near the axis and reaches the speed of light on a cylin- 
drical surface. There is no field in the space outside this 
surface. The field viewed in the co- moving frame of refer- 
ence obeys the force- free relation in the non-relativistic sense 
(V X B) X B = as there is only magnetic field present, 
whereas in the frame of reference at rest with respect to the 
axis the net electromagnetic force by taking into account 
both the interaction of the magnetic fields with the currents 
and the interaction of the electric fields with the charges is 
zero. In the limit of low velocities the equation obtained is 
similar to the Grad-Shafranov partial differential equation 
studied for the case of magnetic fields in equilibrium in the 
absence of pressure, but for relativistic velocities they differ 
considerably. We proceed on the derivation of the equation 
and then we solve it for various cases. 
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2 THE PROBLEM 

Assume a non-resistive plasma of negligible inertia contain- 
ing a magnetic field B and an electric field E in cylindrical 
geometry {R, (j), z), these fields do not depend on z. The elec- 
tric field is related to the magnetic by the definition of ideal 
MHD E — —V X B, where v is the velocity of the plasma 
normalized to the speed of light. In addition assume that 
the plasma expands with 

v = -R. (1) 

ct ^ ' 

The magnetic field can be described by a fiux function P — 
P{v,(f>) for the Br and B^ components, and leave without 
any constraint the Bz component yet, apart from dBz/dz = 
' 



B = \IP{v,4)) X £ + B,(R,(t),t)z , 



or 

„ 1 dP- 1 dP 2 „ - 

R od> R ov 



(2) 



(3) 



The term arising from the gradient of P will be called copla- 



/ - 1 dP \ 

The magnetic and the electric field must satisfy Maxwell's 
equations. Indeed, a magnetic field of the above form is by 
construction divergence-free. The induction equation is 

Vxi. + if 

c ot 



0, 



(5) 



then by substituting equations Q and (j4| into equation (O 
we find that the equations for the R and </!> components are 
satisfied. In order the z component of equation ([5| to be 
equal to zero, B^ has to satisfy the following partial differ- 
ential equation 

we re-express Bz — g{t)W{v,(j}), substituting this into equa- 
tion ([6| gives 



t— -f2 = 0. 



(7) 



therefore g oc t~^, so B^ can be written as Bz — ^14^(11, 0), 
where W is an arbitrary function of v and iji. 

The two remaining Maxwell's equations are used to 
evaluate the charge (p) and the current density (j), 

V -E 



Att 



1 / „ „ dE\ 
j = — cVxB- — 
47r V dt J 



(8) 



(9) 



As we are looking for fields which are force-free in this gener- 
alized sense we demand that the total force the electric and 
the magnetic fields exert to charge and current densities is 
zero. This force is 



F = pE + 



jxB 



(10) 



Thus, by substituting the current and charge density found 



above into equation (|10p we evaluate all three components 
of the force. From the Fz component we take 

then we multiply equation (|lip by (1 — u^)^/^ and re-arragne 
the terms 



d{[i 



,,213/2 



W) dP 



d(j} 



(1 



u ) I — Zvil 

ov 



■ V 



2x1/2 



w 



dP 



dv 
= 0. 



(12) 



(13) 



Setting /3 = (1 - v^f^^W equation gUl reduces to 

9/3ap _ a/3ap _ ^ 

d(j} dv dv d<j) ' 

which is the Jacobian of /3 and P with respect to v and 4>. 
As it equals zero, (3 = f3{P), thus W ^ {1 ~ u^)-^/^/3(P), 
where /3(P) is an arbitrary function of P. Given that, the z 
component of the magnetic field becomes 



Bz = 



i?2(l 



■ u2)3/2 

Using equation (|14|) in the force equation (|10p the R and 
components of the force are both zero when 



(14) 



2n 



+ v{l 



V 



j0 ■ (15) 



This is the cylindrical version of Prendergast's equation, 
which is the basic equation deduced in this paper and 
our major aim hereafter is to find analytical and semi- 
analytical solutions. It simplifies if we apply the transfor- 
mation u — 1/v 

tl 



,2 ,.d'P dP d^P 



du d(jP 



{v? - I) 



(16) 



This is a convenient form of the force-free equation, which 
we shall use later in this paper to find analytical solutions. 



3 SEPARABLE SOLUTIONS 

In what follows we study equation (|15|) and forms that allow 
separable solutions of the form 



P = V{v)^{(f)). 



(17) 



In order to do so, we investigate forms of {3 which permit 
separation of variables. 13 is any function of P and this gives 
a great freedom on the solutions of the basic equation (|15p . 
which may be highly non-linear if such a form is chosen. At 
first we study the form of the equation in the non-relativistic 
limit, then solutions in the absence of a z component of the 
magnetic field and then solutions of the linearized form of 
equation (|15|l . 



3.1 Non-relativistic limit 

For u <C 1 equation (llSp reduces to 



2SPP_ dP_ dPP_ 

dv'^ dv d<P 



-v^0 , 



(18) 
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Separable solutions are allowed in two cases. In the linear 
case for /3/3' — CqP, equation ((18} becomes 



V d f dV\ 2 2 



1 2 



(19) 



where ci is the constant of separation and it is a real number 
as $ has to be periodic. The solution to the angular part is 



<E> = ca cos{ci(j)) + cb sin(ci(; 



(20) 



The differential equation for V admits Bessel functions of 
order ci for solutions 



V{v) = ccJci (cQu) + cdYci (cow) . 



(21) 



The constants ca, cb, cc and cd can be evaluated subject to 
the boundary conditions. The horizontal shear causing the 
z component of the field is parametrized by the constant cq. 

The other class of separabl e solutions is for V = Vbw~' 
and f3l3' = eg P^+^Z' derived by (Il994h . who studied a 
problem of the same geometry in the static case. Indeed if 
we substitute the forms above in equation (|18p . the terms 
involving v cancel and we take an ordinary differential equa- 
tion for $ 



0. 



(22) 



This equation can be solved numerically, subject to the 
boundary conditions and the parameters cq, / and Vq. It is 
possible to by-pass the numerical solution of equation 
and construct an approximate solution for I c lose to zero as 
the one found bv lLvnden-Bell fc Boilvl (|l994l ). by normaliz- 
ing to unity / = ^/^max we find 



l + l 
I 



lnCOsh(d((?!) - 0max)) 



(23) 



= icosh 'expji-. 



where d^ = i±icg(yo$max)'/' and , 

If / is no more small then we can use the more accurate and 
compi cated form described in appendix A of iLvnden-Belll 
l|2006D . 

Although the mathematical treatment described above 
is the same to Aly's case, the two equations describe dis- 
tinct physical configurations. The fields in our problem are 
not static but expand while being force-free. This growth 
can also be conceived as some flux generated from the base 
and expanding with velocity equal to R/ct . In Aly's case 
the fields are static, however it is possible to achieve ex- 
pansion due to the shearing of the field lines on the surface 
of the cylinder. That problem is studied as a sequence of 
static solutions where the transition from one solution to 
the next one takes place slowly enough so that the displace- 
ment current is negligible. In our case the shearing does not 
change explicitly with time as it is imposed from the begin- 
ning. However it is possible to assume the field line s to be 
sheared slowly, indeed as discussed by iLvnden-Belil (l2006l ) 
MHD systems where the boundary conditions change slowly 
have negligible induction currents and pass through a se- 
quence of equilibrium stages. In our case, if the shearing of 
the field lines takes place slowly, then a series of force-free 
solutions of expanding fields describes this process. Thus 
the shearing of the field lines may take place simultaneously 
with the expansion. 



3.2 Relativistic Solutions 

Here we investigate solutions of the fully relativistic equa- 
tion (|15p without neglecting any term. First we study solu- 
tions in the absence oi a, z component for the magnetic field 
which gives a coplanar magnetic field without any current 
and charge density and then solutions containing all three 
components by assuming a linear form for the right hand 
side of the differential equation. 



3.2.1 Current-free solution 

The simplest non-trivial solution of equation (I15|l exists in 
the absence of a z component of the magnetic field, /3 = 0. 
It is a magnetic field that emerges from the surface of a 
cylinder where no shear is imposed. In this configuration 
there is no current or electric charge. We shall use equation 
H16|) as this form is more convenient, let P = V'(u)$(0) and 
derivatives to be denoted by dash, equation (|16|) becomes 



(m" - 1)V" + uV _ _ 2 

The solution of the angular part is of sinusoidal form 
^ — CA sin(ci(^) + Cb cos{ci<l)) , 
the solution for V{u) is 



V = cc{u + y/u'^ - l)"i + coin - y/u^ - 1)" 



(24) 



(25) 



(26) 



This corresponds to a coplanar magnetic field, where every 
magnetic field line lies on a plane normal to the axis of the 
cylinder. 



3.2.2 Linear solution 

Unlike the non-relativistic equation, which admits two 
classes of separable solutions when /? 7^ 0, in the relativistic 
regime, separable solutions exist only in the linear case, and 
there is no analogue to the non-relativistic powerlaws. To 
show this assume separable solutions of the standard form; 
in division of equation (I16|) by P it is 



{u^ - 1)V" + uV ^ 



1313' 



(27) 



the first term of the left hand side is only a function of u 
and can be expressed as F{V), the second term similarly is 
G($), as for the right hand side let H{P) = Pl3'/P. The next 
step is to act with the operators $d/d^ and in the resulting 
expression with the operator Vd/dV; the final equation is 



V 



d{- 



dH 



dV dlnP 



Vu2 ~lj 



Then name Y ■ 



d(lnP)2 ■ 
^ and substitute above. 



d'-H 



dH 



VdY 

YdV ^ ~ d(lnP)2' dlnP 



n . 



(28) 



(29) 



From equation (f29)) it is F oc V" which leads to 1^ = 
Vb( )^'^" and the form of H{P) that permits separation 
of variables is H = ci + C2P~". When substituting these 
forms in the right hand side of equation (|16p it becomes a 
sum of a term that only depends on u and one that only 
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depends on however when the left hand side of the equa- 
tion is considered the terms involving u do not sum up to 
a constant thus it is impossible to separate the equation by 
using power laws. It is only the linear form for /3/3' = CqP, 
which permits separation of variables, indeed 



1)V"+UV' , 2 
r Cg 



V 



(^2 



1 



2 



(30) 



The angular part admits again sinusoidal solutions, as for 
the V{u) the equation to solve is 



(u^ - 1)V" + uV' 2 

V + 



2 
Cl , 



(31) 



which cannot be integrated analytically and we solve it nu- 
merically. 

The numerical solution of equation (|31|l depends on the 
choice of the parameters and the boundary conditions. As- 
sume that some flux emerges from a cylindrical surface Rq 
at time to which in velocity space lies «o = Ro/{cto), the 
boundary conditions for V on this surface are determined 
by the fields. The fields are expressed in terms of V and $, 
they are 



Bf 



R 



R 



B. = 



i?2(l _ ^2)3/2 



(32) 



(33) 



(34) 



The value of V on the surface the flux emerges from is re- 
lated to the intensity of the R component of the magnetic 
field and the derivative V'{vq) to the intensity of 4> compo- 
nent, in what follows we use the normalised value V{vo) — 1. 
Then the parameters cq and Ci are chosen and we integrate 
the differential equation. The physical meaning of the pa- 
rameter Co is related to the horizontal shear of the field lines 
and therefore since Bz is proportional to co, the more shear 
is induced on the base of the field the larger this parameter 
is. An inspection of the differential equation even without 
proceeding to the numerical solution demonstrates the decel- 
erating effect of this term, greater shear leads to a stronger 
Bz component of the magnetic field. This component is pro- 
portional to 7^^*^ and exerts a strong force to the magnetic 
field. Thus, the flux decreases rapidly for i; close to unity. 
After this point the solution undergoes oscillations (Figures 
1, 2) and the field forms disconnected appendages. These 
structures are the result of the linear assumption made pre- 
viously. It is now possible to draw the field lines for any 
time ti (Figures 3, 4) as we have the expressions for V and 
"1> and the field lines are given by substituting expressions 
(l32l). (I33ll and ^ in 



(35) 



dR 



Rd(j) 

Bj, 



dz 



The other parameter ci defines the order of the multipole 
chosen, for instance ci = 1 is a cylindrical dipole and ci = 2 
gives a quadrupole etc. If a non-integer value is chosen, then 
the problem can be studied within a confined arcade, of any 
opening angle, however the solution will not have a physical 
significance outside this domain. Finally the value of the 
derivative V'{vo) is chosen and the differential equation is 



V(u) 




Figure 1. The numerical solutions of equation H30| l. V is nor- 
malized to unity at vq = 0.1, cq = 1 and ci = 1, the solid line 
is for V'{vo) = -10.1, the dashed for V'{vo) = -10, the dot- 
ted for V'{vo) = -9.9 and the dash-dotted for V'{vo) = -9.8. 
This decrease in the absolute value of the derivative leads to a 
weaker at vq, thus the field emerges closer to the perpendic- 
ular and the flux function reaches a greater Vmax before starting 
the oscillations. 



integrated numerically. In the solutions plotted we choose 
four values for it, we find that for a smaller absolute value 
of V'{vo) the flux reaches greater distances in v space; thus 
fields where the ratio Bb.(vq) / B^iva) is bigger have fluxes 
that reach greater distances, so the more radial the flux is 
at the base the further it gets before it starts oscillating. 



4 PHYSICAL QUNATITIES 

In this section we study the shear and the energy stored in 
this conflguration subject to cq. We assume that the same 
radial flux emerges from the surface by keeping V{vo) the 
same and it reaches a maximum velocity i^max also constant, 
where ^(umax) = 0, given these boundary conditions we 
solve equation (|31|l for various values of the parameter cq. 
Then we use the solutions to study the shear and the energy. 



4.1 Shear 

The amount of shear for a given fleld line can be evaluated by 
integrating equation psp after substituting the flelds from 
equations (|32)l34p . An element of a field line labeled by Pi 
that lies between v and v -\- dv m velocity space is sheared 
by 



dz 



coPivdv 



V{1 - (P,/V)2)l/2(1 _„2)3/2 



(36) 



An important feature of this quantity is that it does not 
change as the field expands. We evaluate the total shear 
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Figure 2. The numerical solutions of equation H30| l. V is nor- 
malized to unity at = 0.1, cq = 2 and ci = 1, the solid line 
is for V'{va) = -10.1, the dashed for V'{vo) = -10, the dotted 
for V'{v(,) = -9.9 and the dash-dotted for V'{vo) = -9.8. This 
decrease in the absolute value of the derivative leads to a weaker 
at vq, thus the field emerges closer to the perpendicular and 
the flux function reaches a greater Vraax before starting the oscil- 
lations. Compared to the lines plotted in Figure 1 the z field is 
stronger here and the oscillations start earlier. 




Figure 3. The magnetic field lines. They emerge from the yellow 
cylindrical conductor and they expand radially. The ones reaching 
greater distance from the cylinder have more shear. 




Figure 4. Surfaces of constant P, the line of sight is parallel to 
the axis of the cylinder. The inner circle corresponds to a cylinder 
of radius v = 0.1 and the outer circle to one of radius v = 1. 

of a field line by integrating equation (|36|) from vo where 
the field line emerges, up to the maximum vt it reaches, 
this Vt is given by solving equation V{vt) = Pi, then by 
symmetry the shear from the top to the other foot point is 
the same. Inspection of equation (|36|l shows that the shear 
depends on cq as it appears directly and indirectly in the 
integral through the form of V which is a numerical solution 
depending on the parameter co, there is also a dependence 
on the choice of the field line which also appears directly 
in the equation but also indirectly through the end point 
of integration Vt which depends on Pi. In Table 1 we use 
the results of the linear solution to evaluate the shear {Z) 
of the field line which has been stretched the most for a 
given Cq. The numerical integration shows actually that as 
Co increases the shear increases, a proportionality relation 
between cq and Z is a first approximation, and the deviation 
of this approximation is due to the indirect dependence on 
Co through the other quantities appearing in (|36|l . 

4.2 Energy 

4-2.1 Energy flow 

The total energy carried by an electromagnetic field in a 
volume Vi is 

Etot = ^ /(^' + ■ 

We can express the forms of the electric and the magnetic 
field in terms only of t and v and study the energy contained 
in a volume in velocity space and how it evolves with time. 
By substituting the forms of the magnetic fields from equa- 
tions (|32I34[) and by the fact that E = —v x B the energy 
contained in a volume in velocity space with boundaries vo 



6 K.N. Gourgouliatos 



CO 


Z 


Enp 


En^ 


Etot 


0.0 


0.000 


1.000 


0.000 


1.000 


0.1 


0.175 


1.000 


0.001 


1.001 


0.2 


0.354 


1.001 


0.002 


1.003 


0.3 


0.532 


1.003 


0.004 


1.007 


0.4 


0.738 


1.006 


0.008 


1.014 


0.5 


0.934 


1.011 


0.014 


1.025 


0.6 


1.160 


1.019 


0.024 


1.043 


0.7 


1.426 


1.034 


0.039 


1.073 


0.8 


1.690 


1.064 


0.068 


1.132 


0.9 


2.057 


1.131 


0.131 


1.262 


1.0 


2.332 


1.336 


0.312 


1.648 



5 APPLICATIONS 



Table 1. The physical quantities of the system for cq ranging 
from to 1, evaluated by using the solution of equation II31I I 
subject to the boundary conditions V(O.l) = 1, 1^(0.95) = and 
ci = 1. As Co increases the shearing of the field lines (Z) increases. 
The energy carried by the field also increases, where Eup is the 
energy carried by the coPlanar component of the magnetic field 
and the E^ of the electric; En^ is the energy carried by and 



and Umax, z from to zi and (j) from to 27r, is 

t,tot 



8tv 







(ct)4(l -u2)3 



I V 

^vdvdcj) . 



(38) 



The form of the integral in equation (|38p suggests that the 
energy within a constant volume in velocity space changes 
with time. It is constant only when there is no z component 
of the magnetic field. The Foynting vector of an electromag- 
netic field is 



S = —E xB, 

and by substituting in terms of V and $ we find 



(39) 



S = 



cv f 
4^1 



{ctr 

v{ctY 



,2^3 



R 



^;(ct)3(l-t)2)3/2 



(40) 



As the 7? component of the Foynting vector is positive, the 
energy fiows outwards as the fields expand. 



4-2.2 Energy as a function of the shear 



We have integrated equation (|38p from i; = 0.1 to Wmax = 
0.95 for a range of values of cq at = 1. We used the 
form of V found by solving the differential equation (|3ip 
subject to the boundary conditions stated in this section; 
the results appear in Table 1. The overall conclusion is that 
the energy increases with cq. This is due to two reasons, the 
z component of the magnetic field increases, thus it carries 
more energy. This increase on Bz has also a side effect, for 
the field to remain force-free the coplanar component of the 
magnetic field Bp and = — w x Bp have to increase to 
balance the extra force. Thus the energy carried by them 
becomes larger as well. It is evident from Table 1 that both 



Observations of 7-ray fiares l| Palmer et al.l |2005| ) suggest 
that they are associated to strong magnetic fields expanding 
with relativistic velocities. These flares are thou ght as po- 
tenti al origins of short duration 7-ray bursts C Hurlev et al.l 
l2005h . The relativistic solution of the equations studied in 
this paper can be applied to the initial sta ges of magnetar 
giant flares emerging. It has been proposed l|Lvutikovll2006l ) 
that giant flares from 7-ray repeaters are formed through 
processes similar to coronal mass ejections on the Sun. As 
opposed to solar arcades, they reach high Lorentz factors. 
Our model describes the electromagnetic field of arcades 
which expand and reach relativistic velocities. When the 
flare expands with a great velocity it is essential to take 
into account the relativistic effects as extra terms appear on 
equation p5|) compared to the non-relativistic form, equa- 
tion (|18p . In our study the fields expand only in one dimen- 
sion, this is a reasonable assumption when the magnetic ar- 
cades are not large compared to the radius of the object 
where they emerge from or when there is a primary direc- 
tion of expansion. An other issue is the neglect of inertia 
forces, this is fine when the plasma is underdense and most 
energy is carried by the magnetic field. However very close 
to the speed of hght, no matter how small the density is, 
it will be no more negligible as it will be multiplied by a 
large Lorentz factor. When the expansion in a second dimen- 
sion becomes important then the properties of the spherical 
geometry have to be taken i nto account as they are pre- 
sented bv lFrendergastI l|2005h : [Gourgouliatos fc Lvnden- Belli 
(|2008t ). 



6 CONCLUSIONS 

In this paper w e derived the force- free analogue of Frender- 
gast's equation (|Frendergastll2005l ) in the case of cylidrically 
expanding magnetic fields. These solutions offer a theoreti- 
cal insight to both relativistic and non-relativistic expanding 
systems. The non-relativistic solutions have structural sim- 
ilarities to static ones where time is not explicitly included, 
but it is implied via a series of static solutions where the ex- 
pansion is merely due to shearing of the magnetic field lines 
whereas the fiux emerging from the surface is constant. In 
our case the field expands due to two reasons: shearing of 
the field lines as in the static problem and due to an increase 
of the fiux emerging from the surface. The increase of flux 
appears directly in the equations. 

The plane parallel geometry also sets constraints to the 
applications of these solutions. Astrophysical systems are 
more likely to occur in spherical geometry, a cas e stud- 
ied previously (iGourgouliat os fc Lvnden-Belll |2008| ) , how- 
ever this particular geometry leads to results of satisfactory 
accuracy compared to those of the spherical geometry when 
the radius the magnetic flux extends is comparable to the 
radius of curvature of the surface the magnetic field emerges. 
Systems where the expansion is mainly in one dimension are 
also described by this geometry. An astrophysical system of 
interest is that of solar arcades, this model describes the 
combined effect of the extra flux emerging from the base 
and of the shearing of the field lines provided that the size 
of the arcade is small compared to the radius of the Sun. A 



Relativistically expanding cylindrical electromagnetic fields 7 



detailed comparison of plane parallel mode ls and spherical 
ones can be found in iGourgouliatos I l|2008l ) and the loss of 
equilibrium in these models le ading to the ope ning of the 
field lines has been studied by lUzdenskil (|2002i ). It is also 
shown that the energy contained in the magnetic arcade in- 
creases when more shear is imposed. 

In relativistic systems the demand of simultaneous ex- 
pansion and shearing has some problems as the field at the 
top expands very fast whereas the shearing of the field lines 
at the base has to take place slowly. Thus the message for 
the Bz component to increase will arrive much later at the 
top. For this reason we consider that the magnetic field is 
sheared before the expansion takes place. 
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